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We demonstrate that the form and location (not the size or 
spacing) of the energetically preferred geometrical structure 
of the crystalline quark-hadron mixed phase in a neutron star 
is very sensitive to hnite size terms beyond the surface term. 
We consider two independent approaches of including further 
finite size terms, namely the multiple reflection expansion of 
the bag model and an effective surface tension description. 
Thus care should be taken in any model requiring detailed 
knowledge of these crystalline structures. 
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I. INTRODUCTION 

The possible existence of quark matter in the dense 
cores of neutron stars was proposed in pioneering work 
two decades ago |Q-|^ and has been further studied by a 
number of authors H-p. The theory has recently been 
reexamined and new insights into the nature of the phase 
transition were gained which drastically affect the struc¬ 
ture of neutron stars 0 - In previous work the mixed or 
coexistence phase was strictly excluded from the star and 
a large density discontinuity occured at the radial point 
corresponding to the (constant) pressure of the mixed 
phase. It is now understood that the common pressure 
and all properties of the two phases in equilibrium vary 
as their proportion. Consequently the mixed phase can 
occupy a large radial extent of some kilometers dimen¬ 
sion between the two pure phases. Moreover, the mixed 
phase region is expected to achieve an energy minimum 
by formation of a crystalline lattice of the rarer phase 
immersed in the dominant one. The dimensions and geo¬ 
metrical forms of the lattice will vary with the proportion 
of the phases and hence with depth in the star ||^. The 
discontinuity in density of the early work is replaced by a 
region made of an intricate pattern of crystalline forms. 

The assumption underlying the previous picture is that 
the confining transition from quark matter to hadrons 
is a first order transition with effectively one conserved 
charge, namely the baryon number. This situation arose 
in some papers through treatment of the star as be¬ 
ing made solely of neutrons (which is beta unstable and 
therefore an excited state of the star) and in others by 
the imposition of local charge neutrality on a beta stable 


system which effectively also reduced the phase transition 
to one with a single conserved charge 0 - 

In either case therefore, the transition was made to 
resemble that of the liquid-vapor transition in water, a 
constant pressure-temperature one. However first order 
phase transitions in substances of two or more indepen¬ 
dent components such as neutron stars in beta equilib¬ 
rium (for which baryon number and electric charge num¬ 
ber are the conserved quantities) behave quite differently 
[0 . This follows since the conservation laws are not local 
ones but are global. In fact imposition of local neutrality 
violates Gibbs conditions for phase equilibrium. When 
the two phases are in equilibrium with each other they 
can rearrange the concentration of the conserved charges 
to optimize the total energy while conserving the charges 
overall. In general, energy will be minimized by a differ¬ 
ent concentration of conserved charges at each proportion 
of the phases; hence all properties of the phases also vary 
with proportion. 

In nuclear systems it is the isospin restoring interaction 
(responsible for the valley of beta stability in nuclei) that 
exploits the degree of freedom available by the possibility 
of rearranging concentrations while conserving charges. 
The consequence is that regions of confined and decon- 
fined matter will have opposite charge densities and a 
structured Coulomb lattice minimizes the energy. Else¬ 
where it was found that the unstructured phase would be 
preferred if the surface tension where too large |jll[. We 
have misgivings about this conclusion however, which we 
outline next. 

In Ref. 0 the surface tension was treated as a free 
parameter, which it is not. It was therefore possi¬ 
ble to find a value that was large enough to put the 
structured phase at a higher energy than the unstruc¬ 
tured one. However, opening a degree of freedom has 
the effect of lowering energy, or leaving it unchanged, 
not raising it. Were it possible do an exact calcula¬ 
tion of phase equilibrium including the possibility of 
spatial structure, without introducing the bulk approx¬ 
imation corrected for Coulomb energy, surface energy, 
curvature, etc, the above physical principle would be 
obeyed. It is therefore incorrect in the approximation 
scheme (bulk-|-surface-|-coulomb-|-curvature) to arbitrar¬ 
ily choose a value of surface tension that places the struc¬ 
tured phase above the unstructured one. 
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The form, size, spacing and location in the star of the 
crystalline mixed phase was first calculated by Glenden- 
ning and Pei , where the sum of Coulomb and surface 
energies was minimized as a function of the proportion 
of the two phases in equilibrium. Here we extend their 
calculation to include further hnite size effects, i.e. cur¬ 
vature energy and even higher order terms, for the three 
simple geometries, spheres, rods and slabs. We do this in 
two completely independent ways, and show that the in¬ 
clusion of additional finite size effects has a signihcantly 
larger influence on the preferred geometrical form, than 
on the actual size of the form. 

In Sec. II we give a general description of the mixed 
phase. Sec. Ill and IV are devoted to a description of 
the two different approaches to include further hnite size 
corrections, and Sec. V contains a general discussion and 
a summary of our results. 


II. THE MIXED PHASE 

The equation of state of the conhned hadronic phase 
is calculated as in Ref. HQ- The coupling constants 
of the theory are determined by the bulk nuclear prop¬ 
erties, binding energy, BjA = —16.3 MeV, saturation 
density, po = 0.153 fm“^, symmetry energy coefficient, 
asym = 32.5 MeV, compression modulus K = 240 MeV 
and effective nucleon mass at saturation, m* jm = 0.78. 
The ratio of hyperon coupling to mesons as compared to 
nucleon couplings is chosen in accord with Ref. [Q to be 
Xa = 0.6, X:_j = 0.658. 

The bulk properties of the quark phase is described 
in the bag model at zero temperature with the strong 
coupling constant, as = 0. The relevant expressions for 
the pressure, energy density, baryon number and charge 
density can be found in Ref. Throughout, we use a 
bag constant, = 180 MeV. We take for the quark 
masses niu = 5 MeV, rud = 10 MeV and rus = 150 MeV. 
No heavier quarks are present. While the importance on 
the bulk properties from using non-zero u and d quark 
masses is negligible, the surface tension is increased about 
15 percent in the so-called thin wall limit (only massive 
quarks contribute to the surface tension). 

The relation between the pressure and energy den¬ 
sity, and the neutron star’s Schwarzschild radius is found 
from solving the Oppenheimer-Volkoff equations. Mov¬ 
ing inwards through the star, the mixed phase becomes 
energetically favorable when the energy density has in¬ 
creased to a value where the bulk pressures of the two 
phases have become equal. In our approximation (vol- 
ume-|-surface-|-Coulomb-|-....) the formation of geometri¬ 
cal structures in the mixed phase only occurs if the sum 
of finite size and Coulomb energy densities is less than 
about 10 MeV/fm^. In the next section, where finite size 
corrections are determined from the bag model this does 
not pose a problem as the sum of energy densities is at 


most 7 MeV/fm^. In Sec. IV, where the surface energy 
cannot be explicitly calculated from the model, we ensure 
that the sum does not exceed 10 MeV/fm^ by appropri¬ 
ate adjustment of the surface tension. The geometrical 
structure of the mixed phase occurs to a good approx¬ 
imation against the background of the bulk structure, 
except in the outermost regions of mixed phase, where 
the surface pressure from the geometrical structures are 
comparable with the bulk pressures. The Debye screen¬ 
ing lengths were estimated in Refs. 00 to be about 5 
fm. This is approximately the typical radius of the geo¬ 
metrical structures and thus screening is also only a small 
effect. We will neglect these minor complications. If the 
neutron star is massive enough it will have a pure core 
of quark matter. We focus on a neutron star at the mass 
limit, which with our choice for the equation of state is 
at 1.454M0 and a radius of 10.32 km. 

The mixed phase is subdivided into Wigner-Seitz cells, 
that have total charge zero. Adjacent cells therefore do 
not interact. Each cell has the same form as the one 
structure of rare phase it contains. The surrounding re¬ 
gion of the other phase is treated as a bulk system. For 
the simple discrete geometric forms described in []l6| , |T^ , 
drops, rods and slabs, each characterized by the dimen¬ 
sionality d equal to 3, 2, and 1, respectively, the Coulomb 
energy per volume for the three structures can be written 
as 


Ec = 2 tt [phix) - Pqix)? r'^xfdix) = C{x)r'^, (1) 


where ph and pq are the charge densities of the hadronic 
and quark matter at volume proportion y = VqjV of 
quark matter. The ratio of structure volume to cell vol¬ 
ume is denoted by a:: = {rjTZY. When quark matter is 
the rare phase, a; = y, and otherwise, a: = 1 — y. In the 
case of drops or rods, r is their radius and 77. the half 
distance between centers, whereas for slabs, r is the half 
thickness. The function fd{x) is in all three cases given 

by 


fd{x) 


1 

d + 2 


^(2-dx^-^/'^) 

d-2\ J 


+ X 


( 2 ) 


The preferred geometrical form is found by minimizing 
the sum of finite size and Coulomb energies. We take two 
different approaches to get expressions for the finite size 
corrections beyond the surface term. 


HI. EINITE SIZE CORRECTIONS IN THE BAG 
MODEL 

From the multiple reflection expansion method p^ ] 
used with MIT bag boundary conditions [Q, the finite 
size corrections to the energy can be explicitly calculated. 
This has been done up to curvature corrections for mass¬ 
less quarks, but the ansatz provided by Madsen [|^ for 
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the curvature term for massive quarks has shown to re¬ 
produce exact mode-filling calculations for quarks at zero 
temperature. The use of MIT bag boundary conditions 
imply that the thickness of the boundary wall is zero, the 
so-called thin wall approximation. In this approximation 
the two phases in equilibrium are treated independently, 
i.e. any finite size energy contributions from the hadronic 
phase just add to the quark matter contributions. The 
surface energy per volume for the quark matter is 


Es = 


dxa 

r 


S{x) 


( 3 ) 


and the curvature energy per volume 


_^d{d-l)xj _^r(x) 

— =*=- o -~ ^— ^5 


( 4 ) 


where “-I-” is used when x = x, and ” when x = 1 — x- 
The curvature term vanishes of course for slabs. Expres¬ 
sions for the surface tension, cr, and curvature coefficient, 
7 , can be found in Ref. |^. a decreases from about 37 to 
25 MeV/fm^, and 7 decreases from 25 to 12 MeV/fm as 
the chemical potentials drop out through the star. The 
nuclear liquid drop model give a surface tension ahad — 1 
MeV/fm^ for ordinary nuclear matter in vacuum that, by 
comparison with the quark contribution, can be ignored 
in the thin wall approximation |p 0 | . 

When the quark phase is the rare one, the sum of finite 
size and Coulomb energy densities diverges to plus infin¬ 
ity as a function of r in both the r —> 0 and r —> 00 lim¬ 
its, and thereby ensuring the existence of a local energy 
minimum. But the sign change in the curvature energy 
for a bulk quark phase has devastating consequences, as 
the sum of energy densities diverges to minus infinity for 
vanishing r. First we cannot be sure there is a local mini¬ 
mum, and even worse the energy can be lowered by creat¬ 
ing more and more smaller and smaller hadronic bubbles 
in the quark matter. The sign change in the curvature 
energy is not a pathological defect of the model. The 
curvature energy is proportional to the mean curvature 
integrated over the surface of the geometry, and is there¬ 
fore a signed quantity unlike the surface area. This has 
explicitly been confirmed by Mardor and Svetitsky in 
the high temperature, zero chemical potential limit. We 
take this as a warning of how important it is to ensure 
that enough terms in the expansion are included. 

The next term in a finite size expansion of the en¬ 
ergy, at least for massless quarks, is of the form Zo/r 
for spheres, and is identified as the Casimir or zero-point 
energy. Unfortunately, Zq has not been evaluated in the 
multiple reflection expansion of the bag model. However, 
similar procedures, give Zo ^ 1 1^,^. Bag model fits 
to the hadronic spectrum are considerably improved if a 
zero-point term of — 1.84/r is added to the energy [ p4[ . 
About half of it can be explained in terms of center of 
mass corrections from localization of the hadron . We 


therefore choose to treat Zg as a phenomenological pa¬ 
rameter of order one for all three geometries. For the 
bulk quark phase, where -Ep < 0 , we will insist upon Zo 
being positive; whereas in the bulk hadronic case, where 
Er > 0, we will allow both signs for Zg. Even if Zo is 
negative, the sum of finite size and Coulomb energies will 
not be negative except for radii so small that we do not 
believe in the expansion anyway. This always ensures the 
existence of a local energy minimum. 

For cylinders and parallel plates the zero-point terms 
are Eoyt = Zollr"^ and E^iate = ZoPja?’, respec¬ 

tively. I is the length of the structure, r the cylinder 
radius, and a the distance between the plates, in our no¬ 
tation 2r. The zero-point energies per volume are then 


Ez 


Z{x) ZoX 

j,4 


1 

16’ 

d = 


d = 

TT ’ 


_3_ 

477 ’ 

d = 


( 5 ) 


The sum of finite size and Coulomb energies can be writ¬ 
ten as 


E = C{xy + 


Sjx) ^ Z{x) 


( 6 ) 


The minimization procedure is carried out numerically 
for the three geometries in consideration, under the con¬ 
straint that 2Ec = Es + 2Er + 4:Ez. 



FIG. 1. Mixed phase region in a neutron star showing 
the diameter D (solid curves) and spacing S (dashed curves) 
as a function of the radial Schwarzschild coordinate. The 
curves on top are for massive u and d quarks, in the bottom 
ones they are assumed massless, “h drops” denotes hadronic 
drops immersed in bulk quark matter, and so forth. 


We plot, as a function of the radial Schwarzschild coor¬ 
dinate, R, the diameters of the different geometric struc¬ 
tures present in the mixed phase region of the neutron 
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star, and the spacing between adjacent structures im¬ 
mersed in the bulk phase. The core of the neutron star 
up to i? = 2.0 km consist of pure quark matter, while 
beyond R = 7.9 km the neutron star is composed of or¬ 
dinary hadronic matter. In Fig. ^we have neglected the 
curvature and zero-point energies and show the increase 
in structure size of about 5%, when the u and d quarks 
are not assumed massless. All six geometric forms are 
present in this plot. A new form becomes energetically 
favorable at each non-differential point on the diameter 
curve. The discontinuities are only present because of 
the discrete geometries we treat. Their location is in¬ 
dependent of the surface tension, only the radii of the 
geometries are affected. The radii lie between 2.5 and 
7.5 fm, and therefore screening effects are only of minor 
importance even for the largest structures. 

Fig. ^ is a plot similar to Fig. |] without the spac¬ 
ing curves. Included is now the curvature energy and 
three zero-point energies, assumed positive in both bulk 
phases. In comparison with Fig. two striking new fea¬ 
tures emerge. First, the solid diameter curve, that rep¬ 
resent the smallest zero-point energy, Zq = I, abruptly 
drops at i? ~ 3 km, but without changing the preferred 
geometry, hadronic drops. The reason is that a second 
minimum in the sum of finite size and Coulomb energies, 
Eq. show up at around r = 1 fm, if the zero-point 
term is small enough. 



FIG. 2. Similar to Fig. but without the spacing curves. 
The solid, dotted and dashed curves are for Zo equal to 1, 
1.5 and 2, respectively. Zo is assumed positive in both bulk 
phases. 

For increasing Zg the second minimum will eventually 
not be the global one, this happens when Zg > 1.5, and 
will finally even disappear. Even though the sum of en¬ 
ergies is positive for Zg = 1 as it should be, we neither 
expect nor believe the model to be valid for such small 
radii. So we must take Zg > 1.5. Second, no hadronic 


slabs are present. This is no big surprise, since the cur¬ 
vature energy is zero for the slabs, while negative for 
hadronic drops and rods. Hadronic rods are only pre¬ 
ferred in a very narrow range when Zg is unacceptably 
small. Thus the inclusion of further finite size corrections 
has large implications for the location, and even presence 
of the geometric structures in both bulk phases; whereas 
the sizes are unnoticeably affected, except for small Zg, 
where we do not believe in the model. 

In Fig. which is similar to Fig. |^, we have fixed the 
zero-point energy for the bulk quark phase and show how 
the location and size of the preferred geometry changes 
for three different values of Zg for the bulk hadronic phase 

13- 



R (km) 

FIG. 3. Similar to Fig. |^. The solid, dotted and dashed 
curves are for Zg equal to 1.6, 0 and -1.6 in the bulk hadronic 
phase. Zo = 1.6 in the bulk quark phase. 

For the bulk quark phase the three curves are of course 
identical. Again we see that the structure size is rela¬ 
tively less affected by the change in Zg than the loca¬ 
tion of the geometry. The radial extent of the quark 
drops decrease with decreasing Zg. The division between 
quark slabs and rods does not exhibit the same simple 
behaviour. 

We have assumed that Zg for all three geometries are 
identical, at least in each of the two phases. We have no 
reason to expect this should be correct, but the assump¬ 
tion reduces the possible six different Zg’s to two, still 
enough to show the sensitivity of the location of a geo¬ 
metrical structure to the parameter choice. If the small 
contribution from the hadronic phase is added to the sur¬ 
face energy, only a 3-4% increase in the surface tension, 
and the outcome plotted in Fig. the only noticeably 
difference is that hadronic rods are now preferred in a 
narrow range between hadronic drops and quark slabs, 
the radii are not visibly affected. Furthermore, the small¬ 
est acceptable value for Zg is reduced to 1.3. Fig. shows 
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the reason for this sensitivity. In this figure the sum of 
finite size and Coulomb energies have been plotted as a 
function of the radial coordinate for each of the three 
geometries. The discontinuity at R = 4.6 km is solely 
due to the sign change in the curvature energy. The 
zero-point term is positive in both bulk phases. Notice 
especially how close the rod and drop curves lie in most 
regions. It should be obvious that even small changes in 
parameters could displace the locations of the preferred 
geometries. 



R (km) 

FIG. 4. The sum of finite size and Coulomb energies as 
a fnnction of Schwarzschild coordinate for Zo = 1.6 in both 
bulk phases. The solid, dotted, and dashed curves represent 
slabs, rods, and drops, respectively. 


IV. THE MIXED PHASE IN AN EFFECTIVE 
SURFACE TENSION DESCRIPTION 

In connection with the cosmological phase transition, 
Kajantie et al. have derived an expression for the 
effective surface tension for a diffuse interface between 
quark and hadronic matter in equilibrium, based on a 
simple Taylor expansion of the surface tension of a pla¬ 
nar surface. Also in the high temperature regime the 
curvature term is negative in the bag model for hadronic 
bubbles in bulk quark matter . The sign and magni¬ 
tude has been confirmed by lattice QCD calculations [p^ . 
This inevitably leads to a negative local minimum in the 
free energy even above the phase transition temperature. 
The model of Kajantie et al. has the attracting feature 
of reproducing the negative curvature term when the ex¬ 
pression for the effective surface tension is expanded, but 
still higher order terms prevent the puzzling negative lo¬ 
cal energy minimum from ever forming. 

We cannot immediately use the results for the effec¬ 
tive surface tension presented in that paper. The reason 


being that in our case the pressures and the surface ten¬ 
sion are functions of two variables, namely the baryon 
chemical potential, and the electron chemical poten¬ 
tial, /ie, instead of being just a function of temperature. 
This reflects the presence of two conserved charges in 
our case. But because the electron distribution is uni¬ 
form and /ie <C /in, the effective surface tension can still 
be written in the form, generalized to also include cylin¬ 
drical and planar geometries (see Appendix A), 


cr(r) 


cr(oo) 

l±{d-l)S/r' 


{d-l)6 

r 


< 1 , 


(7) 


where “-I-” is used when x = 1 — x and ” when x = x- 
d is the dimensionality, (t(c») the surface tension of a 
planar surface, and h is a parameter of dimension length 
given as 

( 5 = ( 8 ) 

where nn,z is the baryon number density for each of the 
two phases. 

We take, as in Ref. El , for the surface tension 

cr(x) = const X [e,(x) - eh{x)] A, (9) 



R (km) 

FIG. 5. Similar to Fig. The curves are for S equal to 
0 (solid), 0.5 (dotted), 0.9 (dashed), and 2 fm (dash-dotted), 
respectively. 


where is the energy density for each of the two 
phases. L is a length scale of the order of the interface 
thickness, about 1 fm. The constant is chosen so that 
the sum of finite size and Coulomb energy densities does 
not exceed 10 MeV/fm^. With this choice for the surface 
tension (u varies between 51 and 66 MeV/fm^), S can ac¬ 
tually be calculated, but for our purpose of demonstrat¬ 
ing the large sensitivity of the location of a geometrical 
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structure to the parameter choice, we will treat (5 as a free 
parameter. Estimates based on Eq. show that (5(x) 
varies between a fraction of one fermi to a few fermi’s, but 
in some regions around x — 0-5 it may actually become 
negative. 

We plot in Fig. ^ the structure diameter as a function of 
the radial coordinate for four different values of 5, where 
Eq. is valid. The solid curve where i5 = 0, i.e. no 
finite size terms beyond the surface term, is similar to the 
solid curves in Fig. |^. The maximum structure radii are 
somewhat larger than in the bag model, because of the 
larger surface tension. When 6 is increased, the hadronic 
drops will take up an increasing fraction of the bulk quark 
mixed phase region at the expense of hadronic rods and 
slabs, both of which finally even disappears within the 
range of the chosen i5’s. For the bulk hadronic matter it 
is the quark slabs that take up an increasing fraction at 
the expense of quark rods and drops. If 5 is negative, 
the general tendencies become inverted, quark drops and 
hadronic slabs will eventually dominate the mixed phase 
region. 

Also in this model the structure sizes are relatively 
less affected, than the location and even presence of a 
geometric structure. This is expected because the gen¬ 
eral trends are of course identical in the two approaches 
of including further finite size corrections, as the leading 
contribution comes from the curvature term, which are 
similar in the two approaches. The curvature term cal¬ 
culated from the bag model is reproduced for 5 w 0.6, 
but of course additional terms are different in the two 
approaches. 

V. CONCLUSION 

We have used two different models to include finite 
size corrections beyond the surface term in the sum of 
finite size and Coulomb energies. Within the multiple 
reflection expansion of the MIT bag model the surface 
tension and curvature coefficient can be explicitly calcu¬ 
lated. Because of the negative curvature term, present 
when hadronic matter is immersed in bulk quark mat¬ 
ter, we were forced to add an additional term, which was 
identified as a zero-point energy stemming at least partly 
from Casimir effects. The zero-point term was treated as 
a phenomenological parameter, since terms beyond the 
curvature term have not been determined from the ex¬ 
pansion method. 

The effective surface tension description was originally 
developed for the cosmological quark-hadron transition, 
where a negative curvature term is also present. In the 
present paper this approach was generalized to the zero 
temperature, high baryon and electron chemical poten¬ 
tial limit. For both models the conclusions are qualita¬ 
tively the same. It is the form and location, that relative 


to the size and spacing, become mostly affected by in¬ 
clusion of additional finite size terms. The reason is, as 
shown in Fig. that the minimized sum of finite size and 
Coulomb energies differ very little between the different 
geometries, especially between rod and drop like struc¬ 
tures. Therefore even small changes in the finite size 
terms may significantly displace the location and even 
presence of a given geometry. 

This large sensitivity make calculations involving a de¬ 
tailed understanding of the crystalline mixed phase un¬ 
reliable at present, even for the three simple geometries. 
But it also tells us that the wide range of glitch phe¬ 
nomena observed in pulsars may be closely related to 
restructuring in the solid crystalline region of the pulsar. 
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APPENDIX A: THE EFFECTIVE SURFACE 
TENSION 

In contrast to the cosmological transition studied in 
Ref. the pressures of the two phases and the surface 
tension are now functions of two variables, namely the 
baryon chemical potential and the electron chemical 
potential /Tg. 

By /rjj, /Tg and Pc we denote the values of the chemical 
potentials and pressure when the two phases are in equi¬ 
librium. In the proximity of the transition, the hadron 
and quark pressures can be written as 

Phi^^ni — Pc ^e,/i(/^e Me) (-^f) 

and 

— Pc Me); (-^^) 

where ) is the baryon number density 

in each of the two phases, and riei = {) the 

corresponding electron number density. 

The pressure difference becomes 

Ah = ijln^q nn,/i)(Mn M^) {^c,q ne^/i)(^g Me); 

(A3) 

where the last term vanishes because of the homogeneous 
electron distribution. 
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The surface tension is expanded as 

a{fJ.n,^^e) ^ - Hn) - ^le), 

(A4) 


where a[ is the partial derivative with respect to fii. Since 
Me "C /in we may ignore the third term compared to the 
second term. 

Combining Eqs. ( |A3| ) and (^) with the Laplace con¬ 
dition for mechanical (meta)stability, 


n n (d- l)cr(/i„,/ie) 

— 


(A5) 


generalized to also include cylindrical and planar geome¬ 
tries give a form for the surface tension, Eqs. © and (||), 
similar to the one in Ref. 
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